Fluid squeeze-out between rough surfaces: comparison of theory with experiment 
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We study the time dependency of the (average) interfacial separation between an elastic solid 
with a flat surface and a rigid solid with a randomly rough surface, squeezed together in a fluid. We 
use an analytical theory describing the fluid flow factors, based on the Persson contact mechanics 
theory and the Bruggeman effective medium theory, to calculate the removal of the fluid from the 
contacting interface of the two solids. In order to test this approach, we have performed simple 
squeeze-out experiments. The experimental results are compared to the theory predictions. 



1. Introduction 

Contact mechanics between solid surfaces is the ba- 
sis for understanding many tribology processes [H-ItJ such 
as friction, adhesion, wear and sealing. The two most 
important properties in contact mechanics are the area 
of real contact and the interfacial separation between 
the solid surfaces. For non-adhesive contact and small 
squeezing pressure, the average interfacial separation de- 
pends logarithmically and the (pro jected) contact 
area linearly on the squeezing pressure 111). Here we 
study how the (average) interfacial separation depends 
on time when two elastic solids with rough surfaces are 
squeezed together in a fluid. 

The influence of surface roughness on fluid flow at the 
interface between solids in stationary or sliding contact is 
a topic of great importance both in nature and technol- 
ogy. Technological applications include leakage of seals, 
mixed lubrication, and removal of water from the tire- 
road footprint. In nature fluid removal (squeeze-out) is 
important for adhesion and grip between the tree frog or 
gecko adhesive toe pads and the countersurface during 
rain, and for cell adhesion. 

Almost all surfaces in nature and most surfaces of 
interest in tribology have roughness on many different 
length scales, sometimes extending from atomic distances 
(~ 1 nm) to the macroscopic size of the system which 
could be of order ~ 1 cm. Often the roughness is fractal- 
like so that when a small region is magnified (in general 
with different magnification in the parallel and orthog- 
onal directions) it "looks the same" as the unmagnified 
surface. 

Most objects produced in engineering have some par- 
ticular macroscopic shape characterized by a radius of 
curvature (which may vary over the surface of the solid) 
e.g., the radius R of a cylinder in an engine. In this 
case the surface may appear perfectly smooth to the 
naked eye but at short enough length scale, in general 
much smaller than R, the surface will exhibit strong ir- 
regularities (surface roughness). The surface roughness 
power spectrum C (q) of such a surface exhibits a roll-off 
wavelength Ao << R (related to the roll-off wavevector 
qa = 2tt/Xq) and therefore it appears smooth (except 



for the macroscopic curvature R) on length scales much 
longer than Ao . In this case, when studying the fluid flow 
between two macroscopic solids, one may replace the mi- 
croscopic equations of fluid dynamics with effective equa- 
tions describing the average fluid flow on length scales 
much larger than Ao, and which can be used to study, 
e.g., the lubrication of the cylinder in an engine. This 
approach of eliminating or integrating out short length 
scale degrees of freedom to obtain effective equations of 
motion, which describe the long distance (or slow) be- 
havior is a very general and powerful concept often used 
in physics. 

In the context of fluid flow at the interface between 
closely spaced solids with surface roughness, Patir and 
Cheng [12|, |l3| have showed how the Navier- Stokes equa- 
tions of fluid dynamics can be reduced to effective equa- 
tions of motion involving locally averaged fluid pressures 
and flow velocities. In the effective equation occurs so 
called flow factors, which are functions of the locally av- 
eraged interfacial surface separation. They showed how 
the flow factors can be determined by solving numerically 
the fluid flow in small rectangular units with linear size 
of order (or larger than) the roll-off wavelength Ao intro- 
duced above. In Ref. [14| one of us has developed an 
analytical theory for the pressure flow factors based on 
the Persson contact mechanics model and the Bruggeman 
effective medium theory to take into account the topog- 
raphy disorder resulting from the random roughness. We 
will use this theory in the calculations presented below. 

This paper is organized as follows: In Sec. 2 we briefly 
review the basic equations of fluid dynamics and describe 
some simplifications which are valid for the present case. 
In Sec. 3 we describe the experimental method we have 
used to study the interfacial separation and in Sec. 4 
we compare the experimental results to the theory pre- 
diction. The summary and conclusions are presented in 
Sec. 5. 

2. Theory 

2.1 Fluid flow between solids with random sur- 
face roughness 

Consider two elastic solids with randomly rough sur- 
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FIG. 1: An elastic solid with a rough surface in contact with 
a rigid solid with a fiat surface. 



faces. Even if the solids are squeezed in contact, be- 
cause of the surface roughness there will in general be 
non-contact regions at the interface and, if the squeezing 
force is not too large, there will exist non-contact chan- 
nels from one side to the other side of the nominal con- 
tact region. We consider now fluid flow at the interface 
between the solids. We assume that the fluid is Newto- 
nian and that the fluid velocity field v(x, t) satisfies the 
Navier-Stokes equation: 
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where v = rj/p is the kinetic viscosity and p the mass den- 
sity. For simplicity we will also assume an incompressible 
fluid so that 

V-v = 

We assume that the non-linear term v • Vv can be ne- 
glected (which corresponds to small inertia and small 
Reynolds number) , which is usually the case in fluid flow 
between narrowly spaced solid walls. For simplicity we 
assume the lower solid to be rigid with a flat surface, 
while the upper solid is elastic with a rough surface. In- 
troduce a coordinate system xyz with the xy-plane in the 
surface of the lower solid and the z-axis pointing towards 
the upper solid, see Fig. [TJ The upper solid moves with 
the velocity vo parallel to the lower solid. Let u(x, y, t) be 
the separation between the solid walls and assume that 
the slope |Vu| << 1. We also assume that u/L « 1, 
where L is the linear size of the nominal contact region. 
In this case one expect that the fluid velocity varies slowly 
with the coordinates x and y as compared to the varia- 
tion in the orthogonal direction z. Assuming a slow time 
dependence, the Navier Stokes equations reduces to 



d 2 v 

az A 



Here and in what follows v = (v x ,v y ), x = (x,y) and 
V = (d x ,d y ) are two-dimensional vectors. Note that 
v z » and that p(x) is independent of z to a good ap- 
proximation. The solution to the equations above can be 



v « — z(z - u(x))Vp + -^rvo 
2rj u(x) 



(1) 



so that v = on the solid wall z = and v = Vo for 
2 = u(x). Integrating over z (from z — to z = u(x)) 



gives the fluid flow vector 



- T A_Zvp + - u( x)v 



Mass conservation demand that 
du(x, t) 
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where the interfacial separation u(x, t) is the volume of 
fluid per unit area. In this last equation we have allowed 
for a slow time dependence of u(x, t) as would be the 
case, e.g., during fluid squeeze-out from the interfacial 
region between two solids. 

2.2 Viscosity of confined fluids 

It is well known that the viscosity of fluids at high pres- 
sures may be many orders of magnitude larger than at 
low pressures. Using the theory of activated processes, 
and assuming that a local molecular rearrangement in 
a fluid results in a local volume expansion, one expect 
an exponential dependence on the hydrostatic pressure 
V = r lo ex P(p/Po)j where typically (for hydrocarbons or 
polymer fluids) p n « 10 8 Pa. Here we are interested 
in (wetting) fluids confined between the surfaces of elas- 
tically soft solids, e.g., rubber or gelatin. In this case 
the pressure at the interface is usually at most of or- 
der the Young's modulus, which (for rubber) is less than 
10 7 Pa. Thus, in most cases involving elastically soft 
materials, the viscosity can be considered as indepen- 
dent of the local pressure. In the applications below 
the nominal pressure is only of order ~ 10 4 Pa and the 
pressure in the area of real contact of order ~ 10 6 Pa, 
so that the dependence of the (bulk) viscosity on the 
pressure can be neglected. Nevertheless, it has been ob- 
served experimentally [IB, 16 1, and also found in Molecu- 
lar Dynamics (MD) simulations frR flj| . that the effective 
viscosity r] (defined by a = rjv/u, where a is the shear 
stress, u the separation between the surfaces and v the 
relative velocity) of very thin (nanometer thickness) fluid 
films confined between solid walls at low pressure may be 
strongly enhanced, and to exhibit non-Newtonian prop- 
erties. In addition, for nanometer wall-wall separations, 
a finite normal stress is necessary for squeeze-out, i.e., the 
"fluid" now behaves as a soft solid and the squeeze-out 
occurs in a quantized way by removing one monolayer 
after another with increasing normal stress [liij. In the 
application below we only study the (average) separa- 
tion between the walls with micrometer resolution, and 
in this case the strong increase in the viscosity for very 
short wall separations becomes irrelevant. 



2.3 Roughness on many length scales: effective 
equations of fluid flow 

Equations (2) and (3) describe the fluid flow at the 
interface between contacting solids with rough surfaces. 
The surface roughness can be eliminated or integrated 
out using the Renormalization Group (RG) procedure. 
In this procedure one eliminate or integrate out the sur- 
face roughness components in steps and obtain a set of 
RG flow equations describing how the effective fluid equa- 
tion evolves as more and more of the surface roughness 
components are eliminated. One can show that after 
eliminating all the surface roughness components, the 
fluid current [given by (2)] takes the form 



J = A(u)Vp + B(u)v 



(4) 



where A and B are 2x2 matrices, and where it(x, t) 
and p(x, t) now are locally averaged quantities. In gen- 
eral, A and B depend also on Vp (see Ref. [2CJ]). but 
for the low pressures (and pressure gradients) prevailing 
in the application presented below, we can neglect this 
effect. If the sliding velocity v = and if the surface 
roughness has isotropic statistical properties, then A is 
proportional to the unit matrix and is usually written as 
A = —u 3 (j)p(u)/(12r]). In this case from (3) and (4) we 
get 



du 



zrVp(w), 
12r? 



Vp = 



(5) 



If u(x, t) is independent of x then (5) takes the form 

du _ u 3 (j> p {u) 2 
dt 12r/ P 



(6) 



In Fig. [5] we show (j) p {u) calculated using the Persson 
contact mechanics and the Bruggeman effective medium 
theory fli}. The figure shows the dependence of 4> p (u) on 
the separation u for the two (copper) surfaces used in the 
study below. The green curve shows the large it-behavior 
predicted by Tripp [2l|: 

where (h 2 ) = h 2 . ms . See also Ref. [3] for the calculation 
of higher order corrections. The small contact pressure 
involved in the experiments reported on below results in 
relative large (average) separation between the surfaces, 
u > 1.4/i rms on both surfaces (as calculated using the 
theory developed in Ref. (9L lioj). 

2.4 Fluid squeeze-out 

Let us squeeze a cylindrical rubber block (height d and 
radius R) against a substrate in a fluid. Assume that we 
can neglect the macroscopic deformations of the rubber 
block in response to the (macroscopically) non-uniform 
fluid pressure. In this case u(x, t) will only depend on 




FIG. 2: The fluid pressure flow factor as a function of the 
average interfacial separation u divided by the root-mean- 
square roughness amplitude /i rms . For copper surfaces 1 and 
2. T he g reen curve shows the large u-behavior predicted by 
Tripp [2f|. 



time t and (6) will hold. Equation (6) imply that the 
fluid pressure 



P = Zpti 



u id 



r 

X -R 2 



(7) 



where r — |x| denote the distance from the cylinder axis, 
and where we have assumed that the external pressure 
vanish. pfl U id denotes the average fluid pressure in the 
nominal contact region. Substituting (7) in (6) gives 



du ^ 2u 3 4> p {u)p fiuid (t) 
dt ~ 3r]R 2 



(8) 



If po is the applied pressure acting on the top surface 
of the cylinder block, we have 



PRuid(t) = PO - Pcont(^), 



(9) 



where p C ont is the asperity contact pressure. We first 
assume that the pressure po is so small that for all times 
u » /i rms and in this case 4> P (u) « I. For u >> h rnls 
we also haveQ 



Pcont ~ /3E*exp 



u 
u 



(10) 



where E* = Ej (I — v 2 ) (where E is the Young modulus 
and v the Poisson ratio), and where uq — forms /a. The 
parameters a and 8 depends on the fractal properties of 
the rough surface L 9|. Using (10) and (9) we get from (8): 



dp c 



dt 



2u 3 (p cont (t)) 
3ijR 2 uo 



Pcont (PO -Pcont) , 



(11) 



For long times p C ont ~ Po and we can approximate (11) 
with 



dp c 



_ 2u 3 ( Pa ) 

dt 3rjR 2 UQ 



Po(PO -Pcont) 



Integrating this equation gives 



iWt(*) «P0- [PO - Pcont(O)] exp 



^rms 



where 



(12) 




2hf m j> 2ah 2 msPo ' 



Using (12) this gives 

I , Pcont(O) 

u ~ Moo + I 1 uoexp 

Po 



where Uoo = Ualog(/3E* /po). Thus, u(t) will approach 
the equilibrium separation Uoc in an exponential way, 
and we can define the squeeze-out time as the time to 
reach, say, l.Olwoo, which typically will be a few times 
t 1 = [h lms /u(po)] 3 T. In the experiments performed below 
i] = 100 Pas, R « 1 cm and po « 10 4 Pa. Using that 
uq ps h rms and that for the rough surfaces used below 
/i rms w 50 /im, and u(po) ~ 1.4/i rms we get t' ~ 1000 s. 
Thus we expect squeeze-out to occur in about 1 hour, in 
good agreement with the experimental data (see below). 
For flat surfaces, within continuum mechanics, the film 
thickness approaches zero as f -> oo as « ~ T 1 / 2 . Thus 
in this case there is no natural or characteristic time- 
scale, and it is not possible to define a meaningful fluid 
squeeze-out time. 

At high enough squeezing pressures and after long 
enough time, the interfacial separation will be smaller 
than /i rms , so that the asymptotic relation (10) will no 
longer hold. In this case the relation p CO nt(u) can be 
calculated using the equations given in Ref. [l|. Substi- 
tuting (9) in (8) and measuring pressure in units of po, 
separation in units of /i rms and time in units of r one 
obtain 



du 
~dt 



-a 1 (j> p (u)u 3 (l -Pcont), 



(13) 



where a = h ms /uo- This equation together with the re- 
lation Pcont(w) constitute two equations for two unknown 
(it and Pcont) which can be easily solved by numerical 
integration. 

2.5 Rubber block under vertical loading 

Consider a cylindrical rubber block (radius R, height 
d) squeezed between two flat surfaces. If both surfaces 
are lubricated (no friction) the stress at the interfaces 
will be constant p = po = F^/ttR 2 and the change in 
the thickness Ad (assuming linear elasticity) will be de- 
termined by po = EAd/d. However, if the rubber ad- 
here (or is glued) to the upper surface with no-slip the 
situation may be very different 22 , [23| . If d > R the 
stress at the (lower) interface will again be nearly uni- 
form and Ad will be determined by po = E c s Ad/ d where 





pressure 








FIG. 3: A cylindrical rubber block (height d and radius 7?) 
squeezed against a lubricated substrate (no friction). If d > 
R the pressure distribution at the interface will be nearly 
uniform (left) while if d « R (right) the pressure distribution 
will be nearly parabolic. We have assumed that the upper 
surfaces of the rubber cylinders are glued (no slip) to a flat 
rigid disk. 




FIG. 4: The non-uniform hydrodynamic pressure is highest 
at the center of the contact region and will deform the rubber 
block as indicated in the figure. 



E c ff > E but nearly identical to E. In the opposite limit 
of a very thin rubber disk, d « R, the pressure dis- 
tribution at the bottom surface will be nearly parabolic 
p(r) w 2po[l — ( r /R) 2 ] ( see Fig- El and the effective elas- 
ticity E e g >> E. Experiments have shown that when a 
rubber disk is squeezed against a rough surface, even if 
the rubber disk is very thin, the (locally averaged) pres- 
sure distribution at the bottom surface of the rubber disk 



will be nearly constant 24]. This is because the rubber 
is pressed into the ridges on the rough surface under 
vertical loading, and the hydrostatic pressure becomes 
smaller. We also note that while E c s > E determines 
the change in the thickness of the rubber block, the lo- 
cal elastic asperity-induced deformations at the (lower) 
interface will be determined (to a good approximation) 
by the Young's modulus E. 

In Sec. 2.4 we have shown that a flat cylinder sur- 




FIG. 5: Experimental set-up for studies of fluid squeeze-out 
between surfaces of elastic solids. 

face squeezed against a flat substrate in a (Newtonian) 
fluid gives rise to a parabolic fluid pressure distribution. 
This implies that for a very thin (d « R) elastic disk, 
glued to a flat rigid surface, and squeezed against another 
flat surface in a fluid, we expect the bottom surface of 
the elastic disk to remain nearly flat and the assumption 
made in Sec. 2.4 will hold to good accuracy. However, 
if the rubber block is thick enough (d > R) the bottom 
surface of the block will bend inwards as indicated in Fig. 
21 which will slow down the fluid squeeze-out. 

3. Experimental 

We have studied the squeeze-out of fluids between 
solids with rough surfaces as shown in Fig. [5] In the 
experimental set-up a cylindrical silicon rubber block is 
squeezed against a rough counter-surface in the presents 
of a fluid. The rubber block is attached to a dead weight, 
resulting in the driving force = 13.8 N. This force is 
kept constant for all experiments. We have measured the 
downwards movement of the dead weight as a function of 
time using a digital gauge with a relative position resolu- 
tion of 0.5 /im. In order to slow down the whole process, 
we use a very high viscosity silicon oil (Dow Corning 
200 Fluid, viscosity 100 Pas) and a relative low nomi- 
nal squeezing pressure (about 10 4 Pa). In the different 
configurations we either squeeze an elastic silicon rubber 



block, or a rigid glass block, against smooth (glass) or 
rough (copper) surfaces in order to test different aspects 
of the squeeze-out. The rubber blocks have the radius 
R = 1.5 cm and height d = 1 cm, 0.5 cm and 0.3 cm. 
We use a silicone elastomer (PDMS) prepared using a 
two-component kit (Sylgard 184) purchased from Dow 
Corning (Midland, MI) . This kit consists of a base (vinyl- 
terminated polydimethylsiloxane) and a curing agent 
(methylhydrosiloxane.dimethylsiloxane copolymer) with 
a suitable catalyst. From these two components we pre- 
pared a mixture 10:1 (base/cross linker) in weight. The 
mixture was degassed to remove the trapped air induced 
by stirring from the mixing process and then poured into 
cylindrical casts. The bottom of these casts was made 
from glass to obtain smooth surfaces (negligible rough- 
ness). The samples were cured in an oven at 80° C for 
over 12 hours. The rough copper surfaces where pre- 
pared by pressing sandpaper surfaces against flat and 
plastically soft copper surfaces using a hydraulic press. 
Using sandpaper with different grit size, and repeating 
the procedure many times, resulted in (nearly) randomly 
rough surfaces suitable for our experiment. 

The silicon block was placed in the high viscosity fluid 
with some distance to the rough surface. In order to 
avoid kinetic (inertia) effects the initial separation was 
selected to be very small. The nominal force was applied 
by dropping the dead weight with the rubber block at- 
tached to it. The displacement of the dead weight from 
its starting position was measured as a function of time. 

In Fig. [7] we show the power spectrum of the two 
rough copper surfaces 1 and 2 used in our study. The 
area of real contact (at the nominal squeezing pressure 
ss 2 x 10 4 Pa) as a function of the magnification Q is shown 
in Fig. [5] Note that the area of real contact (i.e., the 
contact area at the highest magnification d or wavevec- 
tor = 9oCi) i s rather similar in both cases (equal to 
A = 0. 016A and 0.013A for surfaces 1 and 2, respec- 
tively) in spite of the rather large difference in the rms 
roughness values (h lms = 42 and 88 /mi, respectively). 
This is due to the fact that the rms roughness is domi- 
nated by the longest wavelength roughness components, 
while the area of real contact is strongly influenced by the 
short wavelength roughness components, which are very 
similar on both surfaces (see Fig. [7] for large wavevec- 
tor). The small contact pressure result in relative large 
(average) separation between the surfaces, u ks 1.4/i rms 
on both surfaces (as calculated using the theory devel- 
oped in Ref. (EEfll). 

4. Comparison of theory with experiment 

In Fig. [9] we show the surface separation as a function 
of the logarithm of time when a glass and a PDMS cylin- 
drical block are squeezed against a flat glass substrate 
in a silicon oil. Also shown is the theoretical prediction 
(lower curve). The cylinder is d = 0.5 cm thick and has 
the diameter D = 2R = 3 cm. As expected, the theory 




FIG. 6: Squeeze-out experimental set-up. A cylindrical glass 
or rubber block is squeezed against a substrate with a smooth 
or rough surface in a fluid. The cylindrical body has the 
height d — 0.5 cm thick and the diameter D = 2R — 3 cm. 
The normal load Fn = 13.8 N and the fluid viscosity r\ = 
100 Pas. The vertical displacement u of the upper surface is 
registered as a function of time. 
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FIG. 7: The logarithm of the surface roughness power spec- 
trum as a function of the logarithm of the wavevector for two 
copper surfaces, 1 and 2, with the root-mean-square rough- 
ness 42 (im and 88 fim, respectively. 




FIG. 8: The calculated relative area of contact A/Ao (where 
Ao is the nominal contact area) when a glass and a PDMS 
block is squeezed against two copper surfaces 1 and 2 with 
surface roughness produced as described in the text. 
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FIG. 9: The surface separation as a function of the logarithm 
of time when a glass and a PDMS cylindrical block is squeezed 
against a flat glass substrate in a silicon oil. Also shown 
is the theoretical prediction (lower curve). The cylindrical 
body has the height d — 0.5 cm thick and has the diameter 
D = 2R = 3 cm. The normal load F N = 13.8 N and the fluid 
viscosity r\ — 100 Pas. 
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FIG. 10: The surface separation as a function of the loga- 
rithm of time when PDMS cylindrical blocks with thickness 
10 mm and 3 mm are squeezed against a rough copper sur- 
face 1 (root-mean-square roughness /i rms = 42 /im) in a sil- 
icon oil. Also shown is the theoretical prediction for a flat 
substrate (dashed curve) and for the copper surface 1 (lower 
solid curve). 



result agree almost perfectly with the experimental re- 
sults for the glass cylinder (no fitting parameters), but 
for the rubber block the (average) separation is larger 
and the squeeze-out slower. We attribute this to tem- 
porarily trapped fluid resulting from the upward bend- 
ing (before contact with the substrate) of the bottom 
surface of the rubber block as in Fig. [4] We define the 
"trapped" fluid volume AV as the fluid volume between 
the bottom surface of the block and a flat (mathematical) 
surface in contact with the block at the edge r = R of the 
bottom surface of the block. Using the theory of elastic- 
ity AV = ttR 2 S with 5 = C'Rp/E cfi , where C is a con- 
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FIG. 11: The same as in Fig. HOI but for a more narrow time 
interval. The thin solid line is the calculated squeeze-out when 
the pressure flow factor 4> p (u) = 1 
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FIG. 12: The surface separation as a function of the log- 
arithm of time when a d — 0.5 cm thick PDMS cylindrical 
block is squeezed against a flat glass substrate (lower curve), 
and against the (rough) copper surface 1 (root-mean-square 
roughness h Ilns — 42 fim) in a silicon oil. 
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FIG. 13: The surface separation as a function of the log- 
arithm of time when a d — 0.5 cm thick PDMS cylindrical 
block is squeezed against the rough copper surface 2 (root- 
mean-square roughness htms = 88 (jm) in a silicon oil. Also 
shown is the theoretical prediction (lower curve). 



stant of order unity. For the rubber block with thickness 
d = 0.5 cm we get E cS « E[l + 0.5 x (R/2d) 2 } w 4 MPa 
and 5 « 40 /im, resulting in an increase in the average 
interfacial separation by ~ 40 /xm, which is consistent 
with what we observe. 

Fig. [TOl and [Til show the surface separation as a func- 
tion of the logarithm of time when PDMS cylindrical 
blocks with thickness 1 cm and 0.3 cm are squeezed 
against the rough copper surface 1 (root-mean-square 
roughness h lms — 42 /zm) in a silicon oil. Also shown 
is the theoretical prediction for a flat substrate (dashed 
curve) and for the copper surface (lower solid curve), as- 
suming that the bottom surface of the rubber disk is 
macroscopically flat. Note that the agreement between 
the theory and experiment is much better for the thin- 
ner rubber disk. This is indeed expected since the fluid- 
pressure induced curvature of the bottom surface of the 
rubber is smaller for the thin rubber disk (see Sec. 2.5). 
But even for the thin rubber disk some fluid-pressure in- 
duced bending of the bottom surface of the rubber disk 
is expected, and we believe this is the main origin for the 
slightly slower squeeze-out observed in the experiment 
as compared to the theory prediction. The thin solid 
line Fig. [TT] shows the calculated squeeze-out when the 
pressure flow factor <fi p (u) = 1. In the present case the 
pressure flow factor is close to unity and this explains the 
relative small difference between using <p p = 1 (thin red 
line) and using the calculated 4> p (u) (from Fig. [2]) (thick 
red line). 

In Fig. [12] we show the surface separation as a func- 
tion of the logarithm of time when the d = 0.5 cm thick 
PDMS disk is squeezed against a flat glass substrate 
(lower curve), and against the (rough) copper surface 1 
(root-mean-square roughness h rms = 42 fim) in a sili- 
con oil. Note that before contact with the substrate, the 
fluid-pressure induced bending of the bottom surface of 
the block is the same in both cases (giving overlapping 
curves for t < 300 s). 

In Fig. [I3]we show the surface separation as a function 
of the logarithm of time when the d = 0.5 cm thick PDMS 
cylindrical block is squeezed against the (rough) copper 
surface 2 (root-mean-square roughness h rms — 88 /im) 
in a silicon oil. Also shown is the theoretical prediction 
(lower curve). Again, the bending of the bottom surface 
of the rubber block results in a slower squeeze-out than 
predicted theoretically assuming a (macroscopically) flat 
bottom surface of the rubber block. 

5. Summary and conclusion 

In this paper we have studied the fluid squeeze-out 
from the interface between an elastic block with a flat 
surface and a randomly rough surface of a rigid solid. 
We have calculated the (average) interfacial separation 
as a function of time by considering the fluid flow using 
a contact mechanics theory in combination with thin-film 
hydrodynamics with flow factors (which are functions of 



the (local) interfacial separation) obtained using a re- 
cently developed theory. We have explained the impor- 
tance of the large length-scale elastic deformations on the 
squeeze out. 

The theoretical results have been compared to ex- 
perimental results. The experiment was performed by 
squeezing cylindrical rubber blocks with different height 
d against rough cooper surfaces in the presents of a high 
viscosity fluid (silicone oil). Changing the height d of 
the rubber block, and also performing additional exper- 
iments with flat against flat surfaces, with combinations 
of rigid-rigid and elastic-rigid, we could show the impor- 
tance of both the large length-scale and asperity induced 
elastic deformation on the squeeze-out. In particular, 
large length-scale deformations of the bottom surface of 
the rubber block resulted in (temporary) trapped fluid 
between the elastic solid and the rigid countersurface, 
which drastically slowed-down the squeeze-out. This ef- 
fect is smallest for the thinnest rubber block, for which 
case we find good agreement between the theory (where 
we have neglected the large length-scale deformations of 
the rubber block) and the experiments. Another mecha- 
nism which drastically slows down the squeeze-out occurs 
at much higher nominal pressure (or load) than used in 
the present experiment. This is due to sealed-off fluid 
in the nominal contact region during contaqct forma- 
tion. This effect occurs when the area of real contact 
approaches « QAAo, where the area of real contact per- 
colate resulting in sealed-off regions of fluid, which may 
disappear only extremely slowly, e.g., by diffusion into 
the rubber. This effect was discussed in Ref. |25j and 
seams to be of importance in many applications involv- 
ing high contact pressures, e.g., it may result in a static 
(or start-up) friction force which slowly increases with 
time even after very long time (say one year). 

The squeeze-out of fluids from the interfacial region 
between elastic solids with rough surfaces is very impor- 
tant in many technical applications (e.g. a tires rolling 
on a wet road, wipers and dynamic seals), and the results 
presented in this paper contribute to this important sub- 
ject. 
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